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This paper examined the relationship between remotely sensed anthropogenic heat discharge and energy 
use from residential and commercial buildings across multiple scales in the city of Indianapolis, Indiana, 
USA. The anthropogenic heat discharge was estimated with a remote sensing-based surface energy bal- 
ance model, which was parameterized using land cover, land surface temperature, albedo, and meteoro- 
logical data. The building energy use was estimated using a GIS-based building energy simulation model 
in conjunction with Department of Energy/Energy Information Administration survey data, the Assessor’s 
parcel data, GIS floor areas data, and remote sensing-derived building height data. The spatial patterns of 
anthropogenic heat discharge and energy use from residential and commercial buildings were analyzed 
and compared. Quantitative relationships were evaluated across multiple scales from pixel aggregation to 
census block. The results indicate that anthropogenic heat discharge is consistent with building energy 
use in terms of the spatial pattern, and that building energy use accounts for a significant fraction of 
anthropogenic heat discharge. The research also implies that the relationship between anthropogenic 
heat discharge and building energy use is scale-dependent. The simultaneous estimation of anthropo- 
genic heat discharge and building energy use via two independent methods improves the understanding 
of the surface energy balance in an urban landscape. The anthropogenic heat discharge derived from 
remote sensing and meteorological data may be able to serve as a spatial distribution proxy for spa- 
tially-resolved building energy use, and even for fossil-fuel C0 2 emissions if additional factors are 
considered. 

© 2011 International Society for Photogrammetry and Remote Sensing, Inc. (ISPRS) All rights reserved. 


1. Introduction 

The urban heat island (UHI) phenomenon is formed when 
higher atmospheric and surface temperatures in urbanized areas 
are observed over the surrounding rural areas (Voogt and Oke, 
2003). Based on the surface energy balance theory, UHI is mainly 
caused by the combination of anthropogenic heat discharge due 
to energy consumption, increased impervious surface area, and 
decreased vegetation and water area (Kato and Yamaguchi, 
2005). Human induced energy discharge has an important impact 
on the urban environment in terms of the surface energy balance. 
Quantification of each heat flux in the energy balance, especially 
the human induced anthropogenic heat discharge and its spatial 
pattern, is important to improve the understanding of human im- 
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pacts on the urban environment, a key issue in global environ- 
mental change. 

The methods used to estimate anthropogenic heat discharge can 
be grouped into three major categories: inventory approaches, 
micrometeorologically-based energy budget closure methods, and 
building energy modeling approaches (Sailor, 2011). For example, 
using a inventory approach Lee et al. estimated anthropogenic heat 
emissions in the Gyeong-In region of Korea in 2002 based on the en- 
ergy consumption statistics data (Lee et al., 2009). Kato and Yamag- 
uchi separated the contribution of anthropogenic heat discharge 
and heat radiation due to solar input to the sensible heat flux by 
using an energy balance method (Kato and Yamaguchi, 2005). 
Heiple and Sailor (2008) estimated the hourly energy consumption 
from residential and commercial buildings at 100 m spatial resolu- 
tion in Houston, Texas using the building energy modeling method. 
Each of these methods has its strengths and limitations. Estimation 
of anthropogenic heat emissions using the inventory approach is 
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based on the real energy consumption data at the spatial scale of 
utility service territories (Ichinose et al., 1999; Klysik, 1996; Lee 
et al M 2009; Sailor and Lu, 2004; Taha, 1997). It is difficult to quan- 
tify the spatial distribution of energy consumption at fine spatial 
and temporal scales. Spatial and temporal downscaling can only 
be achieved when additional data such as land use are employed. 
Energy balance method has been applied to estimate energy fluxes 
in urban areas from site and neighborhood scale to city scale (Belan 
et al., 2009; Masson et al., 2002; Oke, 1988; Oke et al., 1999; Pigeon 
et al., 2007; Pigeon et al., 2008). High spatial resolution anthropo- 
genic heat discharge can be estimated based on micrometeorolog- 
ically-based energy budget closure method by the combined use 
of remote sensing and meteorological data (Bastiaanssen et al., 
1998a,b; French et al., 2005; Kato and Yamaguchi, 2005; Kato 
et al.,2008; Schmugge et al., 1998; van derKwastet al., 2009). How- 
ever, as the estimation of anthropogenic heat discharge is based on 
the residual of other components in the surface energy balance, 
each component introduces uncertainties and propagates errors to- 
ward the final estimated anthropogenic heat flux. By integrating the 
geospatial data and simulated temporal pattern of energy con- 
sumption for representative buildings, energy consumption in all 
buildings can be estimated using the building energy modeling 
method (Heiple and Sailor, 2008; Ichinose et al., 1999; Zhou and 
Gurney, 2010). However, the representative buildings may not be 
easily categorized and the use of representative buildings may 
introduce bias in the estimation of spatial and temporal patterns 
of energy use for some buildings as these buildings may have differ- 
ent behaviors of energy consumption in terms of magnitude and 
temporal pattern compared to the representative buildings. 

Although anthropogenic heat discharge was studied using differ- 
ent methods in previous studies, there was limited research for 
cross-examination of the spatial pattern of anthropogenic heat dis- 
charge. Even though previous studies have provided insights into 
the human impacts on the urban environment directly or indirectly, 
simultaneous estimation of anthropogenic heat discharge using 
multiple methods will help to reduce the uncertainties in the study 
of urban heat fluxes and improve the understanding of the role of hu- 
man activities in the surface energy balance and its contribution to 
UHI. This study aims to improve the understanding of the urban sur- 
face energy balance and clarify the spatial pattern of energy related 
human activities with the help of remote sensing data, technologies 
and building energy simulation in a highly urbanized area. Specifi- 
cally, anthropogenic heat discharge and energy use from residential 
and commercial buildings were investigated using two different 
methods in the metropolitan Indianapolis, Indiana region. First, 
anthropogenic heat discharge was estimated using a surface energy 
balance method. Then energy use from residential and commercial 
buildings was calculated using a building energy modeling method. 
Finally, the relationship between remotely sensed anthropogenic 
heat discharge and building energy use was examined across multi- 
ple scales from pixel aggregations to census block groups. 

2. Study area and datasets 

2 A. Study area 

In this study, Indianapolis/Marion County, Indiana was chosen as 
the experiment area (Fig. 1 ). Indianapolis is the capital of the State of 
Indiana, and was listed as the 14th largest city in the USA in 2008, 
with a population of 798,382 (US Census Bureau, 2009). The city 
has a temperate climate without pronounced wet or dry seasons, 
but has obvious seasonal changes. The total annual heating degree 
day value is 5521 F, and the cooling degree day value is 1042F 
(National Weather Service, 2002). The synoptic weather condition 
during the study period is typical of mid-latitude summer time. 
The city of Indianapolis lies on the White River at its confluence with 



Fig. 1 . The study area of Indianapolis City, Marion County, Indiana, USA. 


Fall Creek, near the center of the state. It presently extends into nine 
Marion County Townships, including Pike, Washington, Lawrence, 
Wayne, Center, Warren, Decatur, Perry and Franklin. The city is 
located on a flat plain and is relatively symmetrical, having possibil- 
ities of expansion in all directions, and the city center has expanded 
into adjacent agricultural and rural lands through encroachment 
(Weng and Hu, 2008). Indianapolis is a highly urbanized city with 
spatially heterogeneous human activities, which serves as an ideal 
area for the purpose of this study. 

2.2. Datasets 

Anthropogenic heat discharge at noon time on June 16th, 2001 
(Saturday) was estimated using a surface energy balance modeling 
approach (Fig. 2). The details of the approach were discussed in the 
methodology section. The major data used in the energy balance 
method include meteorological data, land cover, albedo and land 
surface temperature (LST). Impervious surface area (ISA) was also 
developed and utilized as a mask in the analysis. Meteorological 
data including air temperature, shortwave radiation, and wind 
speed at noon time on June 16th, 2001 were retrieved from 3- 
hourly MAT3FXCHM data product at the Modeling and Assimilation 
Data and Information Services Center (http://disc.sci.gsfc.nasa.gov/ 
MDISC/). The data processing of land cover, albedo, LST and ISA is 
discussed below in details. 

Land cover data from US Geological Survey National Land Cover 
Dataset (NLCD) 2001 were collected to retrieve parameters for the 
estimation of surface heat fluxes (Fig. 3a). The land cover classifica- 
tion was achieved by using a classification and decision tree method 
with Landsat imagery and ancillary data (Homer et al., 2004). Main 
land cover types in the study area include open water, urban, ever- 
green forest, deciduous forest, shrub, grassland, cropland and wet- 
land. Urban is the dominant land use in Indianapolis and accounts 
for more than 3/4 of the total area. This dataset was utilized to create 
a look up table between land cover type and modeling parameters, 
including coefficient for ground heat flux, roughness length and 
emissivity, which were developed in previous studies (Kato and 
Yamaguchi, 2005; Kato et al., 2008; Snyder et al., 1998). Land cover 
based parameters used in the estimation of heat fluxes will be fur- 
ther discussed in the methodology section. 
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Fig. 2. Methodology used to estimate anthropogenic heat discharge and energy use from residential and commercial buildings. 
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Fig. 3. Model input data in the study area of Indianapolis including land cover (a), LST (b), albedo (c), and ISA (d). 


The advanced spaceborne thermal emission and reflection radi- 
ometer (ASTER) on-demand L2 surface kinetic temperature data 
with 90 m resolution of June 16th, 2001 were purchased from 


NASA (Fig. 3b). The algorithm for converting ASTER thermal infra- 
red measurements to LSTs has been reported by the ASTER temper- 
ature/emissivity working group (Gillespie et al., 1999), with which 
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surface kinetic temperature was determined by applying Planck’s 
Law using the emissivity values from the temperature-emissivity 
separation (TES) algorithm. LST values calculated using this proce- 
dure are expected to have an absolute accuracy of 1-4 I< and rela- 
tive accuracy of 0.3 K, and surface emissivity values are expected to 
have an absolute accuracy of 0.05-0.1 and relative accuracy of 
0.005 (Gillespie et al., 1999). 

Land surface albedo map in 30 m spatial resolution was gener- 
ated from Landsat surface reflectance and anisotropy information 
based on concurrent moderate resolution imaging spectroradiom- 
eter (MODIS) 500 m observations through the algorithm developed 
by Shuai et al. (2011) (Fig. 3c). The level 1 Landsat 7 enhanced the- 
matic mapper plus (ETM+) image collected on July 11th, 2001 with 
good quality was first converted to 30 m land surface reflectance 
through Landsat ecosystem disturbance adaptive processing sys- 
tem, an automatic calibration and atmospheric correction process 
(Masek et al., 2006). Then the classification-based 500 m Bidirec- 
tional reflectance distribution function (BRDFs) for the homoge- 
nous area were extracted from the concurrent MODIS BRDF 
product (Schaaf et al., 2002) with the highest inversion quality 
(Shuai et al., 2008), followed by the calculation of spectral albe- 
do-to-reflectance ratios for each class at the MODIS pixel scale. 
Next, the spectral Landsat albedos were inversed by applying these 
ratios into Landsat surface reflectance. According to the validation 
with ground measurements, the Landsat albedo meets the accu- 
racy required by regional climate modeling. Compared with the 
MODIS 500 m products, it also provides a detailed landscape tex- 
ture and improved accuracy at the validation stations (Shuai 
et al., 2011). It was assumed that there was no significant change 
of albedo between June 16th and July 11th, 2001. 

ISA was estimated on June 16th, 2001 using the method devel- 
oped by Weng et al. (2009) (Fig. 3d). Through this method, imper- 
vious surface was calculated based on the relationship between 
the reflectance of two endmembers (high- and low-albedo) and 
the reflectance of the impervious areas. By examining the rela- 
tionships between impervious surfaces and four endmembers, it 
is found that impervious surfaces were located on or near the line 
connecting the low- and high-albedo endmembers in the feature 
space. An estimation procedure was thus developed based on this 
relationship by adding the fractions of high- and low-albedo 
endmembers. 

Energy use from commercial and residential buildings was esti- 
mated using a building energy simulation method (right panel of 
Fig. 2). In this method, the simulation of building energy use was 
driven by the typical meteorological year data. Base on data avail- 
ability noon time of June 15th, 2002 (Saturday) was chosen to rep- 
resent building energy use for the typical weekend during the 
summer time to be comparable to the remotely sensed anthropo- 
logic heat discharge. The major data used in the building energy 
modeling includes a GIS database of building footprints, the County 
Assessor’s parcel database, and the building height data. A local GIS 
database of building footprints was collected and used to retrieve 
the floor area of each building. The County Assessor’s parcel data- 
base was a Microsoft Access based dataset. It was used to retrieve 
the building type and age and to link with the GIS building foot- 
print. Building heights were estimated from a remote sensing- 
based digital surface model and a digital elevation model and used 
to calculate floor numbers for each building. 


3. Methodology 

3A. Estimation of anthropogenic heat discharge 

In this study, anthropogenic heat discharge was estimated 
based on a surface energy balance modeling method (Kato and 


Yamaguchi, 2005; Kato et al., 2008; Oke, 1988; Xu et al., 2008). 
The method is illustrated in the left panel of Fig. 2. The surface en- 
ergy balance due to surface properties and anthropogenic heating 
in the near-surface in an urban area can be expressed as follows; 

0/ = H + LE + G - R n (1) 

where Qf is the anthropogenic heat discharge (W m -2 ), H is the sen- 
sible heat (W m _2 ), LE is the latent heat flux (W m -2 ), and G is the 
ground heat flux (W m -2 ), and R n is the net radiation (W m -2 ). 
Based on the surface energy balance model, anthropogenic heat dis- 
charge Qf can be estimated as the residual term. 

In the energy balance equation, net radiation (R n ) is the sum 
of short and long wave radiation at the land surface (Kato and 
Yamaguchi, 2005). It is calculated as 

Rn = (1 - oc)Rs + SaSRd - sR u (2) 

where R s is the shortwave radiation (W m -2 ) which is retrieved 
from the meteorological data, a is the surface albedo which is de- 
rived from the remote sensing data, and the creation of 30 m spatial 
resolution albedo data is discussed in the data section, s is the sur- 
face emissivity, and it is determined based on land type from NLCD 
land cover data (Snyder et al., 1998). s a is the atmospheric emissiv- 
ity, and it is estimated using an empirical equation used by Kato and 
Yamaguchi (2005) 

fe \ V 7 

8 " =i ' 24 6y < 3 > 

where T a is the atmospheric temperature (K) from meteorological 
data, and e a is the atmospheric water vapor pressure (hPa), which 
is estimated using the method by Allen et al. (1998) 

n _ RHmean £°(^min) "T 6 (T max ) ( ^ 

e “"To6 2 (4j 

where e°(T min ) is the saturation vapor pressure at daily minimum 
temperature (hPa), e°(T m ax) is the saturation vapor pressure at daily 
maximum temperature (hPa), and RH mean is the mean relative 
humidity (%). 

R d and R u used to calculate net radiation (R n ) are the downward 
and upward blackbody radiation (W m -2 ), and they are calculated 
using the following equations (Xu et al., 2008) 

R d = gT a (5) 

R u = aTt (6) 

where a is the Stefan-Boltzmann constant (5.67E-08 W m -2 K -4 ), 
and T s is the surface kinetic temperature (K). The process to derive 
the surface temperature was discussed in the data section. 

Sensible heat flux (H) is estimated based on the following equa- 
tion (Allen et al., 1998; Kato and Yamaguchi, 2005; Kato et al., 
2008; Xu et al., 2008) 

T - T 

H = pC„ 5 l-L± (7) 

' a 

where p is the air density (kg m -3 ), C p is the specific heat of air at 
constant pressure (J kg -1 K -1 ), and T 0 is the surface aerodynamic 
temperature. As T 0 is difficult to obtain, remote sensing derived sur- 
face temperature was used to approach T 0 (Kato et al., 2008). r a 
(sm _1 ) is the aerodynamic resistance which is calculated using 
the following simple relation (Allen et al., 1998) 

In in 

r a = V 0m 4 — (8) 

k u z 

where z m is the height of wind measurements (m), z h is the height of 
humidity measurements (m), d is the zero plane displacement 
height (m), z 0m is the roughness length governing momentum 


Y. Zhou et al/lSPRS Journal of Photogrammetry and Remote Sensing 67 (201 1 ) 65-72 


69 


transfer (m), z oh is the roughness length governing transfer of heat 
and vapor (m), k is the von Karman’s constant (0.4), and u z is the 
wind speed (m s _1 ) at height z from the meteorological data. 

Latent heat flux ( LE ) is estimated using the following equation 
(Kato and Yamaguchi, 2005) 


LE = pC P 


e s Ca 
W + r s ) 


( 9 ) 


where e s is the saturation vapor pressure (hPa), y is the psychromet- 
ric constant (hPaK -1 ), and r s is the stomatal resistance (sm _1 ) 
which is estimated based on land type from the NLCD land cover 
data (Kato et al., 2008). 

Ground heat flux (G) is estimated using the following empirical 
equation by Kato and Yamaguchi (2005) 


G — C g R n 


( 10 ) 


where C g is a fixed coefficient based on the NLCD land cover type 
and season (Kato and Yamaguchi, 2005). 

3.2. Estimation of energy use from residential and commercial 
buildings 

Energy use from residential and commercial buildings was esti- 
mated using a building energy modeling method, which is illus- 
trated in the right panel of Fig. 2. The method for C0 2 emissions 
quantification developed by Zhou and Gurney (2010), in which 
on-site fossil-fuel C0 2 emissions from residential and commercial 
buildings were quantified by downscaling the county-level com- 
mercial and residential fossil fuel C0 2 emissions, was modified 
and extended. Using the modified method energy use from resi- 
dential and commercial buildings was estimated for each building. 
First, a building typology was developed to define the building 
types for residential and commercial building stock in Indianapolis. 
A total of 30 building prototypes, 22 commercial and eight residen- 
tial, were defined based on the building type, size, and age with the 
help of the Marion County Assessor’s parcel database. 

Second, hourly building Energy Use Intensity (EUI) - the energy 
used per unit floor area - was simulated using the Quick Energy 
Simulation Tool (eQUEST) for each building prototype (US Depart- 
ment of Energy, 2009). The major parameters used in eQUEST sim- 
ulation including building shell characteristics (e.g., R value), 
internal load (e.g., lighting) and Heating, Ventilating, and Air Con- 
ditioning system characteristics were derived from Huang et al. for 
the North Central US census region (Huang et al., 1991). Indianap- 
olis was specified as the location in eQUEST simulation to retrieve 
the typical weather information for this region. The EUIs and the 
ratio of the non-electric to electric EUI values from the eQUEST 
simulation were calibrated using the US Department of Energy’s 
Energy Information Administration data for the East North Central 
Census Division from its residential energy consumption survey 
and commercial building energy consumption survey (US Energy 
Information Administration, 2001 ; US Energy Information Admin- 
istration, 2003). 

Finally, energy use in each building was estimated from the cal- 
ibrated non-electric and electric EUIs, building area and number of 
floors as follows: 


^building (i) — (EUI tjbuilding ^ + EUI t)building ^ J^building(i)NFbuilding(i) 


( 11 ) 


where EUI® building(i) is the electric energy use intensity of building 
type t for building i, EUI^ e ing(j) is the non-electric energy use inten- 
sity of building type t for building i, A bui i ding(I) is the footprint area of 
building i (which was retrieved from a local GIS database of building 
footprints), and NF bU ji din g (/) is the number of floors for building i 
which was estimated based on building height calculated from a re- 
mote sensing-based digital surface model and a digital elevation 


model using the method by Zhou and Gurney (2010). With the com- 
bination of building area from GIS footprint data, building type from 
the Assessor’s parcel data, floor numbers from building height data, 
and electric and non-electric EUIs for building prototypes, energy 
use in all residential and commercial buildings was estimated in 
the study area. 

4. Results 

4.1. Anthropogenic heat discharge 

The anthropogenic heat discharge in the study area at noon 
time on June 16th, 2001 was estimated using the energy balance 
method. To reduce the uncertainties in the evaluation of anthropo- 
genic heat discharge, especially in less developed areas, ISA data 
were used to exclude pixels with the ISA less than 25%. The result 
of anthropogenic heat discharge estimation is shown in Fig. 4. The 
result shows the spatial pattern and magnitude of anthropogenic 
heat discharge in the City of Indianapolis in a typical summer time. 
The anthropogenic heat discharge ranges from 0 to 400 W m -2 , 
with a mean value of 32 W m -2 . The high anthropogenic heat dis- 
charge occurs in the dense residential and commercial areas such 
as the center of the city and the residential areas in the northwest 
corner of the city. The township of Center has the highest anthro- 
pogenic heat discharge of 78 W m -2 as it has high density of com- 
mercial buildings while the Franklin township with large rural 
areas observes the lowest value. There are negative values in the 
anthropogenic heat discharge estimation, although it theoretically 
should be greater than 0 W m -2 . These negative values are re- 
garded as estimation errors, and were set to 0 W m -2 in the final 
result. 

4.2. Building energy use 

Using the building energy modeling approach, the energy use 
from residential and commercial buildings in each building at noon 
time on June 15th, 2002 was estimated in the study area. The spa- 
tial distribution of the energy use from residential and commercial 
buildings at a re-sampled 90 m spatial resolution is illustrated in 



Fig. 4. The spatial distribution of anthropogenic heat discharge. 
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Fig. 5. To be consistent with the anthropogenic heat discharge esti- 
mation, the building energy use was examined with the exclusion 
of pixels with the ISA less than 25%. The energy use from residential 
and commercial buildings ranges from 0 to 85 W m -2 with a mean 
value of 5.5 W m -2 . The spatial variation of building energy use is 
mainly determined by building density, height and type. The spatial 
pattern shows that the highest building energy use occurs at the city 
center, and that the township of Center possesses the highest energy 
use of 36 W m -2 . The building energy use is also high in dense resi- 
dential areas associated with suburban development in the outer 
ring of Highway 465 in Indianapolis as shown in Fig. 5. 

4.3. Spatial relationship, uncertainty and scale effect 

Although energy use from residential and commercial build- 
ings represents only one component of anthropogenic heat dis- 
charge, it shows an obvious consistency in terms of spatial 
distribution with remote sensing-derived anthropogenic heat 
discharge. Visual comparison between anthropogenic heat dis- 
charge (Fig. 4) and building energy use (Fig. 5) shows that they 
both are higher in the dense residential and commercial areas 
and lower in the rural areas. In order to evaluate the quantita- 
tive relationship two variables were compared at pixel levels 
as well as at the census block level. First, the building energy 
use and anthropogenic heat discharge data were re-sampled to 
a variety of spatial resolutions of 90 m, 180 m, 270 m, 360 m 
and 450 m. Furthermore, the pixels with 0.5% of the ranked 
building energy use in the both ends, which are considered as 
outliers, were excluded in the analysis. The relationship between 
these two variables across these scales and the additional census 
block level is shown in Fig. 6a-f. The result indicates that the 
correlation between building energy use and anthropogenic heat 
discharge is statistically significant at all spatial scales 
(p< 0.001), and that the relationship between two variables var- 
ies with scales. The explained variance of building energy use by 
anthropogenic heat discharge increases with the decreasing spa- 
tial resolution. 

A potential application of this study is to reduce the uncertainty 
in urban energy balance study through independent estimation of 


anthropogenic heat discharge and building energy use. It is worth 
noting that there are uncertainties in both datasets of anthropo- 
genic heat discharge and building energy use which might contrib- 
ute to the discrepancy between these two variables to be discussed 
in the next paragraph. The uncertainty in the estimated anthropo- 
genic heat discharge mainly came from each component in the sur- 
face energy balance. The input parameters of meteorological data, 
land cover, land surface temperature, and albedo may introduce 
the uncertainty towards the final estimation of anthropogenic heat 
discharge. The uncertainty in the estimated building energy use 
mainly came from EUIs using eQUEST simulation and building 
height extracted from a digital surface and a digital elevation mod- 
el. The uncertainty analysis of fossil-fuel C0 2 emissions which is 
based on non-electric building energy use can be found in the 
study by Zhou and Gurney (2010) using a Monte Carlo approach. 
The uncertainty of non-electric energy use in each town is approx- 
imately within 10% of the total energy use. 

The discrepancy between building energy use and anthropogenic 
heat discharge is caused by several factors. First, anthropogenic heat 
discharge includes not only energy use from residential and com- 
mercial buildings, but also other sources such as energy use from 
industrial buildings and transportation. The quantitative relation- 
ship confirms that building energy use accounts for part of anthro- 
pogenic heat discharge. Second, energy use from residential and 
commercial buildings was estimated based on the typical weather 
information at noon time on June 15th, 2002 while the anthropo- 
genic heat discharge was estimated at noon time on June 16th, 
2001. The different dates of the input data may have introduced 
inconsistency between two variables. Third, as only low resolution 
meteorological data such as wind speed and air temperature used 
in the energy balance method, the spatial variation of impact from 
meteorological factors was not considered due to data limitation. Fi- 
nally, image registration of data used in the estimation of two vari- 
ables may also bring in uncertainties in the investigation of the 
relationship. It is one of the reasons that the correlation between 
anthropogenic heat discharge and building energy use varies with 
the spatial resolution. With the decreasing spatial resolution, the er- 
ror of the spatial registration between two images decreases. There- 
fore the correlation between building energy use and anthropogenic 
heat discharge increases as indicated by Fig. 6. 



Fig. 5. The spatial distribution of energy use from commercial and residential 
buildings. 


5. Implications 

There are still challenges to obtaining accurate estimation of 
each component of heat flux in the urban energy balance such as 
the ground heat flux. Accurate estimation of anthropogenic heat 
discharge will help to better understand the urban energy balance 
and separate the contribution of anthropogenic heat discharge 
from that of ground heat flux. By combining with energy use from 
other sources such as industrial buildings and transportation, the 
spatially-resolved energy use from residential and commercial 
buildings will improve the estimation of anthropogenic heat dis- 
charge, and therefore, help to reduce the uncertainties in the study 
of urban energy balance. 

The relationship between building energy use and anthropo- 
genic heat discharge indicates that anthropogenic heat discharge 
may be able to serve as a spatial distribution proxy for the spa- 
tially-resolved energy use in an urban landscape. In addition, be- 
cause of the relationship between energy use and fossil-fuel C0 2 
emissions, anthropogenic heat discharge, which is relatively easy 
to obtain with the help of remotely sensed data, provides the po- 
tential to develop spatially-resolved fossil-fuel C0 2 emissions from 
building sectors with additional considerations. Together with 
point source and mobile emissions, the spatially-resolved fossil- 
fuel C0 2 emissions dataset will help to improve the scientific 
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Fig. 6. The relationship between anthropogenic heat discharge and energy use from residential and commercial buildings at 90 m (a), 180 m (b), 270 m (c), 360 m (d), 450 m 
(e) spatial resolution, and the US census block level (f). 


understanding of the contribution of human induced emissions to 
atmospheric concentration through transport modeling in the ur- 
ban domain. 


6. Conclusions 

In this study, the anthropogenic heat discharge and energy use 
from residential and commercial buildings were estimated using 
two independent methods simultaneously in the City of Indianap- 
olis in the summer time. The anthropogenic heat discharge was 
estimated using a remote sensing method, while the energy use 
from residential and commercial buildings was estimated with a 
GIS-based simulation method. The mean anthropogenic heat dis- 
charge is 32 W m -2 while the energy use from residential and com- 
mercial buildings is 5.5 W m -2 . The visual comparison and 
quantitative analysis show consistency between anthropogenic 
heat discharge and building energy use in terms of the spatial pat- 
tern and magnitude. The spatial relationship between anthropo- 
genic heat discharge and building energy use varies with the 
spatial scales, and the explained variation of the building energy 
use by anthropogenic heat discharge increases with the decreasing 
spatial resolution. The results of the quantitative relationship also 


confirm that the energy use from residential and commercial 
buildings represents only partial anthropogenic heat discharge 
which includes additional components of energy use such as those 
from industrial buildings and transportation. 

The spatial consistency between anthropogenic heat discharge 
and energy use from residential and commercial buildings indi- 
cates that remotely sensed anthropogenic heat discharge may be 
able to serve as a spatial proxy to derive the spatial pattern of en- 
ergy use, and even fossil-fuel C0 2 emissions with consideration of 
additional factors. As the spatial relationship between building en- 
ergy use and anthropogenic heat discharge varies with spatial 
scales, appropriate spatial resolution should be identified in future 
studies to aid the development of spatially-resolved energy use 
datasets. 

The independent estimation of building energy use and anthro- 
pogenic heat discharge will also help to improve the estimation of 
other heat fluxes and the understanding of energy balance in the 
urban areas. In order to reduce the uncertainties in the study of ur- 
ban energy balance, further studies of energy use from other hu- 
man activities, especially industrial activities and transportation, 
should be combined with the estimation of energy use from resi- 
dential and commercial buildings to obtain a more complete pic- 
ture of human energy use in an urban domain. 
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